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1 Introduction 

In this paper we are concerned with a mathematical model leading to a generalised setting of 
relativistic dynamics. 

As it is well-known, within the framework of the special relativity theory, a point particle of a 
constant rest mass mo > moving in an electromagnetic field {A, V) on is usually described 
by the action integral Cdt with the Lagrangian 



(1.1) C = -moc\jl-^^ + -Av-eV, 
c c 
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providing motion equations 



1.2 - " )=-{v^ rot A) - - ^ - e grady, 

at \ L _ v^/ c c at 

with the Lorentz force, T^, on the right-hand side. More generally, it is known (and rather 
surprising) that equations 

(L3) |-(-^)=/ 

are variational "as they stand" if and only if the force on the right-hand side is a Lorentz-type 
force, i.e., = J-l, for some functions A,V on see also [T1])0 Mathematically this 

result means that Lorentz-type force interactions are the only admissible interactions for a massive 
special-relativistic particle, compatible with a variational principle for curves in M.^ parametrised by 
time t £R. On the other hand, there appeared hypotheses that also different kinds of interactions 
should be possible: a significant prediction of such a non-Lorentz-type force is due to Dicke [S], 
however, a formula as well as exact arguments are missing up to now. Hence, a question arises on a 
proper general mathematical setting, providing equations of motion for particles in the Minkowski 
space-time, and clarifying admissible forces. 

A problem of this kind can be tackled by tools of the modern calculus of variations. The above 
variational principle is apparently too restrictive: one should better consider a variational principle 
for one-dimensional submanifolds in M"*. However, there is another interesting possibility, within an 
even more general framework, proposed and studied in |16j : One can utilize a variational principle 
for curves in the Minkowski space-time (M'*, g) together with the relativistic constraint on the four- 
velocity: from the mathematical point of view this is a nonholonomic constraint. The constraint 
defines a genuine evolution space, Q, which is a submanifold of codimension one. 

Then dynamics proceed in Q and are governed by so-called reduced Chetaev equations [501 US] 
that become the desired equations of motion for the problem in question. In the above mentioned 
paper |16j we applied the procedure directly to "conventional" massive particles, moving in an 
electromagnetic and scalar field, and obtained a generalised formula for the force T covering also 
"Dicke-type" interactions. 

In this paper we show that ideas of this kind can be applied to study all kinds of point particles 
admitted by the special relativity theory. We propose a unified mathematical setting suitable not 
only for "conventional" particles (with a positive square of mass) , but also for tachyons (particles 
with a negative square of mass), and particles with zero mass, moving in an electromagnetic and 
scalar field. The idea of [16 to tackle this problem as a variational problem with a nonholonomic 
constraint remains, however, the Lagrangian and the constraint are appropriately modified: the 
Lagrange function we propose is universal for all the kinds of particles (not containing the parti- 
cle's mass), however, particles are distinguished by the constraint (mass appears in the constraint 
condition). In the next two sections of the present paper we explain our approach in general, in 
the rest we study in detail each kind of the particles separately. As the main result we obtain the 
corresponding equations of motion, and formulas for admissible forces. We discuss the meaning of 
the constraint condition, and of the new forces (different from the usual Lorentz force type term), 
arising due to the nonholonomic constraint, and find among them a force complying with Dicke's 
predictions. In particular, we pay attention to equations of motion of tachyons and of massless 
particles, and investigate the corresponding dynamics. In this context one should mention at least 

^Of course, in this general setting, the meaning of functions A, V need not be just a vector and scalar potential of 
an electromagnetic field. It includes also e.g. forces of inertia. 
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two striking properties of massless particles: First, it turns out that dynamics of these particles are 
singular, not obeying the Newton's determinism principle, i.e. these particles belong to mechanical 
systems with internal constraints studied (in the Lagrangian case) by Dirac. Second, we found 
out that, if massless particles are admitted to move in a scalar field, this field does not influence 
on their dynamics. The last section is devoted to a summary of results of this paper, where also 
other interesting properties of massive and massless particles, obtained within our approach, are 
recalled and discussed. 

2 Mathematical background 

From the point of view of mathematics, the main idea is to join variational approach, represented 
by the calculus of variations on fibred manifolds (see 0[TD] or |12j ) with nonholonomic geometric 
mechanics in jet bundles due to [12], [H] (for alternative approaches see e.g. pi 151 151 ITSl 
and others). 

Throughout the paper we shall use Einstein summation convention. Summation over greek 
(respectively, latin) indices proceeds from 1 to 4 (respectively, 1 to 3). 

First, let us recall a necessary mathematical background, adapted to our situation. Consider 
the manifold with canonical coordinates {q"^), 1 < cr < 4, endowed with the Minkowski metric 
field 

(2.1) g = Skidq''®dq^-dq^®dq\ 

Curves in will be represented by their graphs, i.e. sections of the fibred manifold tt : M x M'* — > M. 
If c : s — * c(s) is a curve in K'* parametrised by a real parameter s S M, then the corresponding 
section is the mapping 7 : M 9 s ^ 7(5) = (s,c(s)) = {s,q°'{s)) G M x M''. We shall need the 
manifold 

(2.2) = M X TM-* = M X X 

(the first jet prolongation of tt), called evolution space. We shall consider it with fibred coordinates 
(s, g"^, g"^); note that the "dot" coordinates are defined by 

(2.3) g-o Ji7= A(5-o7), 

ds 

for every section 7 of tt. 

On J^TT we have the so-called contact structure, locally generated by basic contact 1-forms 

(2.4) Lj" = dq'^ - q^ds, 1 < cr < 4. 

In what follows, we shall often consider vector fields and 1-forms, defined on the Minkowski 
space-time R"* (i.e., on the typical fibre of the fibration tt), and call them contravariant and covari- 
ant four-vector fields, respectively. In fibred coordinates a contravariant, respectively, covariant 
vector field reads 

d 

(2.5) u — u'^ — — , respectively, (f> — (f>^dq'^ , 

oq" 

where the components and (^jy , 1 < < 4, are functions on W^, i.e. depend upon the coordinates 
{q"). Due to the structure of the evolution space J^tt (being a product of M and the tangent space 
TW^ = X R"*), contravariant four-vector fields are sections of the bundle TW^ R^, i.e., it 
holds q^ o u = q'^ and 

(2.6) cf ou^v!'{q^,q^,q^,q^), I < v < A. 
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Covariant four- vector fields are sections of the cotangent bundle T*Mf^ M^. 

A Lagrangian is defined to be a differential one-form on the evolution space J^tt, horizontal 
with respect to the projection onto the base M. In fibred coordinates, a Lagrangian is expressed 
by A = L(s, q'^ , q") ds, and action (over an interval [a, b] C M) is a function 



(2.7) S':r(7r)9 7^ ( J^-f*A^ [ {L o J^j) ds e R, 

Jla.b] Ja 



..b] 

where r(7r) denotes the set of sections of tt with domains of definition containing the interval [a, b], 
and J^7 is the first jet prolongation of 7, i.e. a section of J^tt defined by J^j{s) = (s, c(s), dc/ds) = 
{s,q-is),q-is)). 

It is important to notice the geometric meaning of the Euler-Lagrange equations of a La- 
grangian defined on J^tt. Extremals of A are integral sections of a distribution on the evolution 
space J^TT (by a distribution we mean a subbundle of the tangent bundle to the manifold J^tt). This 
distribution, called dynamical distribution, arises as the characteristic distribution of the 2-form 

(2.8) a ^ AaUj" A ds + Ba^uj" A dq" , 

where 

^ ■ -* " ^ ^ ds dq" ^ dq^dq" ^ ' ""^^WW' 
In this notation, Euler-Lagrange equations of A read 

(2.10) + = 0, 1 <a<4, 

or, in intrinsic form, J^'j*i^a = for every vector field ^ on J^tt. The rank of the dynamical 
distribution may be greater than one. It is equal to one if and only if the matrix {B^^,) is regular, 
i.e. if the Lagrangian is regular; in this case trajectories in the evolution space are integral curves 
of a single vector field. For more details on the structure of solutions of Euler-Lagrange equations 
and corresponding integration methods for both regular and non-regular Lagrangians we refer to 
[12] and [E]. 

In this paper we shall consider Lagrangians where the Lagrange function L does not depend 
explicitly upon the parameter s, as relevant from the physical point of view. In this case we have 

(2.11) (L o Ji7)(s) = l{c{s), ^) = L{q%s),r{s)) = {Louo c){s) = (L o u){c{s)), 

hence, we can consider a Lagrange function defined on M x TM^, of the form L = L{u). 

A nonholonomic constraint is a submanifold Q of J^tt, fibred over MxM"^. Whenever convenient, 
we denote by t the canonical embedding of Q into J^tt. We shall consider a nonholonomic con- 
straint of codimension 1. Such a constraint is given by one first-order partial differential equation 

(2.12) /(,s,g^g")-0, 
such that 

(2.13) mnkidf/dq") = 1. 

By the rank condition, we may assume that the constraint is defined by equation "in normal form" 

(2.14) g" = /^(s,9^g'), 
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where 1 < < 4, and 1 < / < 3 (i.e. that f = - h). The given fibred coordinates induce on the 
submanifold Q adapted coordinates {s,q^ .q*^ ,q^). 

A constraint in the evolution space J^tt gives rise to a constraint structure on J^tt, generated 
by the l-form 

(2.15) ^^fds+^u\ 

As shown in [131 [IB], due to the constraint structure, the constraint Q is naturally endowed with 
a distribution (subbundle of the tangent bundle to Q) called canonical distribution, annihilated 
by the one-form 

(2.16) (p = i*ip^ -TTTT^^' + (dg* - /ids) 

Oq^ 

on the manifold Q. The canonical distribution represents admissible directions in the evolution 
space (a nonholonomic "principle of virtual displacements" ) . 

There are two models for describing a constrained system, both will be useful for our further 
considerations. The first one is more traditional, describing the constrained system as a defor- 
mation of the original unconstrained system due to a constraint force, naturally generated by the 
constraint. The deformed system is thus a new mechanical system defined on J^tt. It is represented 
by the 2-form ac = a — ^, where a represents the original mechanical system and 

df 

(2.17) $ = A A ds = A —^Lu" A ds 

is the constraint force, called Chetaev force. Constrained trajectories are then integral curves of 
the characteristic distribution of the 2-form ac, passing in the manifold Q. Equations for these 
curves, called Chetaev equations [3], depend upon one Lagrange multiplier A (to be determined), 
and read 

,r..c.. dL d dL Bf 

2.18 -r = ^ 1< o- < 4. 

^ ' dq'^ dsdq'^ dq'^ ~ ~ 

The second model is more geometrical. It represents a non-holonomic mechanical system as an 
object defined on the constraint Q. Hence, the manifold Q has the meaning of the evolution space 
for the constrained system. Concretely, the constrained system is represented by the class of 2- 
forms [a] — i*a + [0\ on Q, where [if] is the ideal generated by the constraint form if. Now, 
admissible trajectories are integral sections of the canonical distribution, and equations of motion, 
called Chetaev reduced equations are equations for admissible trajectories that are integral curves 
of the characteristic distribution of a. Keeping notations used so far, they take the form 

(2.19) Ai + Buq'^O, l<l<3, 

where 

(2^20) A, = {a, + A,^+ U. + B.,^^- + 



(2.21)A.^(a. + 2B„|L+B«||L)„,. 

It is important to stress that if the constraint is non-integrable (and this will be our case) the 
non-holonomic equations of motion do not come from the "constrained Lagrangian" L ^ L o l = 
L{s,q%q\h). 
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3 A geometric setting for SRT particle dynamics 



In this paper we propose a model for particle dynamics on the Minkowski space, suitable for a 

unified treatment of all kinds of point particles ("classical" - with positive square of mass, tachyons 

- with negative square of mass, and particles with zero rest mass). 

This model treats a particle as a Lagrangian system on a fibred manifold tt : M x M"' ^ M (where 
is endowed with the Minkowski metric), with a non-holonomic constraint in the evolution space 

J^TT — M. X TM^, and is based on the following two axioms: 

(1) The Lagrange function is universal for all particles, and is polynomial in the 
"four- velocity" 



L{u) — — \g{u, u) + (f>{u) — ip — — \gauu"u^ + (pau"' — ip 

(3.1) 

2 I 



where (pa-, 1 < f < 4, and are functions on the space-time M"*. 

Here "universality" means that the Lagrangian does not contain the particle mass. Note that 
coefficients in the Lagrangian have the physical meaning of admissihe fields: behind the (2,0)- 
tensor field g, the Lagrangian contains a covariant four-vector field (p and a scalar field "0. 

(2) A particle is determined by a non-holonomic constraint condition on the "four- 
velocity" : 

(3.2) g{u,u) = 7Wc^ 

where Ai — Ai{q^,q^, q'^ , q^) is a function on the space-time M^. 

The value of the function A4 may vary from point to point; at a point in it may be positive, 
negative or zero. 

Note that if we assume M be continuous then the following property holds: If A4{x) 7^ at a 
point a; G K'' then there is an open neighbourhood U around x such that ^ on [/. This means 
that if M is positive (negative) at a point, it is positive (negative) in a certain neighbourhood of 
that point. 



Since g(u,u) — {{q ) ~ ^p^iiq^) ) ° condition (3.2 I defines a nonholonomic constraint 

(3.3) Q: {q^f -j^m^ = Mc^ 
p=i 



in J^TT = M X TW^ if the rank condition (2.131 is satisfied, i.e., if 

(3.4) rank(-4^-g^-'?^4') = l• 

Hence we have to exclude from M x TM^ the points where m = (the zero section in TM^) . However, 
we need more, namely, the constraint Q be expressible in normal form 

(3.5) q^^h{s,q%q\eA'')- 
This is the case if 

(3.6) gVO. 

Excluding the hyperplane Tio of points in J^i: where = we get a disconnected open suhmanifold 
of the evolution space J^tt, consisting from two connected components, J^tt^ and J^7r_, of points 
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where > and 5* < 0, respectively. On the components we have the induced global fibred 



coordinates (s, q*^, g"^), and the constraint condition (3.2 1 determines a constraint submanifold 
Q C J^n\Ho = Q+ U Q_, defined by the following equations in normal form: 



(3.7) Q+ C JV+ : q^ 



\ 



Mc2 + ^(gP)2 



= 1 



9* = 



A 



Mc2 + ^(gP)2. 

p=i 



(3.8) Q_ C J^7r_ 

In what follows we shall choose for the evolution space of the constrained system the manifol 



Lagrangian (3.11 gives rise to four Euler-Lagrange equations 



(3.9) ^-^^ = 
^ ' dq'' dsdq" 

for graphs of curves c(s) — {q'^{s)) in the space-time M'', having the following form: 



(3.10) 







dtp 


dql 


dq" ) 


dq^ 


d4>a 




dtp 


dq^ 


dq" J 


^ dq^ 



-0, 
= 0. 



Equations we are looking for concern motions proceeding in the evolution space Q+. Substitut- 



ing 



(3.11) 



and 



d(t>a 




dip 


dq'- 


dq'' J 


dq'- ' 


dcpa 


904 \ 


dtp 


dq^ 


dq" J 


^ dq^ 



J, Bij^-Sij, l<l,j<3, 

Bil = Bn = 0, i?44 = 1, 



(3.12) h 



\ 



Mc2+^((jP)2. 

p=l 



into (2.201 and (2.211 we obtain 



Ai = q' 



d(pi \ dt/j 



dq' 9g' / dq' 



d(f)i 904 \ dtp 
9g* 9q* 



(3.13) 



Xc2+^((jP)2 
p=l 



dq' dq''' 



B, 



q'ql 
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The desired equations are three equations for sections of the fibred manifold Q+ over M, i.e ., fo r 
curves c(s) = (g'^(s)) passing in the space-time and satisfying the constraint condition (3.7), 
as follows: 

(3.14) Bijq^ =-Ai, l<l<3. 

Now, we have to distinguish two cases: 

3.1 = at a point x E Q+ 



If A4{x) = then the matrix {Bij) is singular at x, hence motion equations {3.14) cannot be put 
into the normal form; explicitly they read 



q'-q^ 



dq'- dq^ 



q^q'- 



(3.15) 



d(t>i _ 904 
2 \dq* dq' 



904 d<j)i \ dip dip 



dq^ dq^ J dq^ dq^ 



where 1 < ^ < 3. Extremals of the constrained problem are curves 7(5) = (s, q'^{s)) satisfying the 
above motion equations and the equation of the constraint 2+ 



(3.16) q" = 



The dynamics proceed in the evolution space 

We shall deal with these equations in more detail later when we shall study massless particles. 



3.2 A< 7^ at a point x e Q+ 

If M.{x) ^ then the matrix (-B;-,) is regular at x, and, due to continuity, it is regular in a 
neighbourhood U around x. Hence, on U, equations (3. 14-) can be put into the normal form 

(3.17) q^ = P. 

The inverse matrix — {B^') takes the form 



/ Mc^ + iq^ 



(3.18) B-^ 



• 1 -2 

q q 



• 1 • 2 

q q 



q^q^ 



q'q' \ 



hence we obtain 



A, 



(3.19) 



5'^q' 
1 



dq'- dq' 



SI I 



\ 



Mc-^ + J^iq^y 
p=l 



dcpA _ d(h 
dq' dq''' 



7Wc2 + ^(gP): 
p=l 



dip 
dq^ 
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Proposition 3.1. (Equations of motion: ^ 0, four-dimensional observer.) Let g be the 

Minkowski metric, (j) a covariant vector field, and tp a function on M^. Extremals of a Lagrangian 
system defined by the Lagrangian 



(3.20) L = -y{u, u) + cPiu) - ^ = -i ((g4)2 - Y^iq^f) + q''cj>^ - V 

p=i 

on M X TM^, and subject to the constraint 

(3.21) g{u, u) = Mc^, > ^ 0, 

are curves 7(5) = {s,q'^{s)), satisfying the following system, of mixed second and first order differ- 
ential equations: 



q' = q ' 



(3.22) 



dq^ dq'- J \^ 



+ 



Mc"^ + ^(gP)2 



d(f)4 d(t>j 
dq} 9g4 



dip 



dip 



dqi Mc^ y dq^ ^ 
(equations of motion) and 



Mc^ + Y^iqpy 



p=i 



dtp 
d^ 



1 < j < 3 



(3.23) 



A4c2 + ^(gf )2 

p=i 



(equation of the constraint). 

The dynamics proceed in the evolution space Q+ C M x TM^ defined by the above constraint 
equation. 

Choosing appropriate coordinates, we can express the above equations in a form adapted to a 
three-dimensional observer. 
Let us denote 

(3.24) g4 = ci 

(time coordinate), and consider on J^n\Ho C M x x new coordinates (s, q^t, v\ q^), defined 
by the transformation rule 

(3.25) = - v'-q^, 1 < / < 3. 

c 

Note that (s, g', t, ,q^) are global coordinates, however, no longer fibred coordinates for the orig- 
inal fibration tt. The meaning of the new coordinates is the following: {t,q^,v^) are coordinates 
on R X ]R3 X 



adapted to the fibration M x 



of the manifold M^, the fibre of our fibred 



manifold tt : M x ^ M; note that {q^,q^, q'^) are cartesian coordinates on R^. In vector notations 
^= (9^)9^)9^)) and V = {v^,v'^,v^) is the usual velocity. 

In these coordinates the constraint Q is given by equation 



(3.26) {l-^){qy = Mc\ 
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Since we assume ^ 0, we obtain the evolution space Q+ expressed by equation 
(3.27) = 



1 



where V — yJ2p^ii'^^)'^ = yJ2p^ii^^^/^^)^ usual three-dimensional speed. We can also 

see that along every constrained path, 



(3.28) ^ = '-q^ 
as c 



M 



1 



Now, equations of motion (3.221 can be transformed eliminating the parameter s as follows: 
First, we have for j = 1, 2, 3 



(3.29) = -(q^) 
Next, if we denote 



1 d 4 1 dt d .4 



' M d_ 

1-4 



' M 



(3.30) (j) = cjj^dq'' = (j)idq^ + (j)4cdt = -Adf- eVdt, 



I.e. 

(3.31) {^l)l = i^2,3 



V, 



where e is the particle charge, we get from (3.191 the force 
(3.32) f=e 



M 



1 



1^ r 

- V X rot A grad V 

c c at 



V d^ 



2 - t;2 



jrad ^. 



Proposition 3.2. (Equations of motion: M ^ 0, three-dimensional observer.) Let (j) 

(3.301 he a one-form (covariant four-vector field), and ij) a function on the Min kowsk i spac e-tim e 
W^. Consider on (R x TW^)\Hq adapted coordinates {s,q'' ,t,v\q'^) defined by (3.24) and (3.25). 
Extremals of a Lagrangian system defined by the Lagrangian A = Lds, 

„.2 s 



(3.33) L 



1 



1 



(9^ 



1 /e 



c \c 



-Av- eV]ct -il) 



and subject to the constraint 
(3.34) g4 = 



1 - ^ 

are curves c(s) = (i(s), q'(i(s))) in M"*, satisfying the following system of differential equations: 

/l - " 



(3.35) A 
dt 



M 



V , 



1 



1 ^1 dA 

e ( - u X rot A grad V 

. c c at 



M 



jrad %l) 



c2 — dt 



(equations of motion) and 



(3.36) 



dt 



M 



(equation of the constraint). The dynamics proceed in the evolution space C M x TK^ defined 
by equation (3.34). 
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We shall see later that equation of the constramt (3.36) has the meaning of mass equation, or, 
if multiplied by c^, of energy equation. 

In what follows, we shall denote by the Lorentz force, 

/ 1 1 dA 

(3.37) Tl = e(-vxrotA grad V 

\c c at 

3.3 The constraint force 

Chetaev equations of motion contain the constraint force <&. It is defined on J^tt = M x TW^ and 
depends upon one Lagrange multiplier A (to be determined). For 



(3.38) f = q'- 
it reads 

3 

(3.39) <i=-Y. 



\ 



Xc2 + ^(gP)2 



^3 



:dgf' Ads + Adg^Ads. 



Chetaev equations take the form 



(3.40) 



if + cf 







dt/j 


.dq' 


dq" J 










9g4 


dq" J 


^ dqi 



-Ag' 
Mc^ + EUiqPy 



= A. 



Multiplying the j/-th of the above 4 equations by q'^ and summing over v = 1,2, 3,4, we obtain 

dq" dq" 



(3.41) {qY-J2q'q') + rr 



1=1 



= A 



Solutions c(s) = {q'^{s)) of Chetaev equations satisfy the equation of the constraint, and conse- 
quently, also its derivative. 



1 _d 

(3.42) q'^^ 2d. 



4-4 



and then ij q 



1 d 



1=1 



^3 f ;.n\2 2 ds 

Hence, the first term in brackets on the left-hand side of equation (3.41 1 is g^(^-^c2), the second 
term is identically zero because of skewsymmetry of the expression m brackets, and the third term 
equals —dijj/ds. In this way, along solutions of Chetaev equations passing in the evolution space 
we have 



(3.43) ^(iTWc^-V') = A 



/Mc^+j:;^,{qpr 

If 7^ 0, we can compute the multiplier A: 



(3.44) A 



Mc? 
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Writing $ = ^o-d?'^ A ds and substituting the obtained expression for A into formula (3.391 for <I>, 
we get components of the constraint force along solutions of the constraint equations as follows: 



(3.45) 



$4 



Mc^ ds 



in coordinates {Syq"^ ,q'^), and 

d 



ds 



<i>i = <i>i = 



(3.46) 



<I>4 



C$4 



r,2 _ y2 

C2 d 



W2 dt 



in coordinates {s,q\t,v'' ,q'^). 

For the sake of simplicity we introduce the following vector notation, that takes into account 
the formula for the constraint force and the relation between equations (3.22) and (3.35) (note 
that they are connected via a multiplier) 



(3.47) Tc^-\ 



M 



i<;<3 



dt 



and call J-c the induced constraint force. Note that J-'c is defined on the evolution space Q+, and 
vanishes whenever w = or 

(3.48) ?/; — = constant of the motion. 



The constraint force appears in the motion equations (3.35): 



Proposition 3.3. If Ad ^ Q then equations of motion (3.35) can he expressed as follows: 
d 



(3.49) 



dt 



M 



1 



= Tl + Tc- 



1-g 
M 



^rad i}} 



dM 



dt ' 



where Tj^ is the Lorentz force and J-'c is the induced constraint force. 

In the rest of the paper we shall be interested in the case Ai — const . We shall discuss 
separately the cases > 0, < 0, and M — 0. Note for constant A4 the above equation takes 
the form 



(3.50) 



dt 



M 



1 ^ 



1-g 
M 



^radi/). 



4 Particles with positive square of mass 

Let us suppose > and constant on M*. We set 
(4.1) M ^ ml. 

As above, we consider either fibred coordinates {s,q'^,q"'), 1 < ct < 4, or adapted coordinates 
{s,q\t,v\q^), 1 < / < 3, on (M X TR^Hq. 
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w 




Figure 1: Projection of the evolution space for 7W > 0. 



4.1 The evolution space 

The constraint Q C (M x TR'^)\Ho is given by equation 

(4.2) iqr-j2i<i''r-rnlc', resp. (q^ffl-^-l)^rnlc\ 

p=l \ c / 

The latter implies — > 0, i.e. 

(4.3) V <c. 

Thus, the speed of a relativistic particle with A4 > is lower than the light 
For the evolution space 2+ we then have 

(4.4) q'= 



1 - % 

where mg > 0. We can see that the constant mp has the meaning of the rest mass of a particle. 
We also denote, as usual, mass and (kinetic) energy by 

(4.5) m — — , c = mc . 
1 - 4 

We can "visualise" the constraint Q and the evolution space in the "velocity space" as 
follows: For illustration we supress one dimension by considering q'^ = and denote {q^, q^ , g^) = 
c{x,y,w). Then Q is given by 

(4.6) 4-4-4=1 

ttiq TOq mg 
(Fig. 1). 
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4.2 Four-momentum 

Since M"* is equipped with the Minkowski metric field g which is a regular symmetric (2, 0)-tensor 
field on M^, every contravariant four- vector field ii on is canonically associated with a covariant 
four-vector field (a 1-form) g{u, ■) on M^. We put 

(4.7) p = g{u,-) 

and call p the four-momentum associated with u. In fibred coordinates, where the section u has 
components vl^ — cf o ii, 1 < a < A, we have p — p^dq'^ , where 

(4.8) Pa ^ Qauu" = g„i,{(f ou), i.e. {pcy)i<cy<i = {-q^ .-cf' .q^) ° u. 

For simplicity of notations, we shall denote components oi p associated with u simply by 

(4.9) {Mi<.<i = {-q\-q\-q\q^). 

where the g's are components of ii. 

In adapted coordinates (s, 9', i, w', q^) we similarly have p — pidq^ + pidt, 

(4.10) p, = -^, pi = cq\ 

c 

Note that, in the latter case, the corresponding "contravariant components" are 

(4.11) = p^ = cq\ 

c 

We denote p= (?j')i<i<3' ^^^^ write 

(4.12) (p.)i<.<4 = (- — , cg4), (p-)i<,<4 = ('^, cg^). 

~ ~ c ^ ^ c 

With help of four-momentum, equation of the constraint Q, i.e. g{u,u) = to§c^, reads 

(4.13) p^p'^ = mgC^, or p,p' -f ^^4^" = mpC^. 
On the evolution space Q+, where 

(4.14) q — — -j^^^= = mc = — , 



components of the four-momentum take the form 

(4.15) {pa)i<cr<i^ i-mv, £/c), {p'')i<a<4 = {mv, £/c), 



Finally, with help of (4.141 we can compute the energy; we obtain the familiar formula 



(4.16) £ = c Y riT^c^ — pip'- = c Y rriQC^ + p^ , 

where p^ = p ■ p= X]f=i(?'')^ denotes the usual square of length of the three-vector p. 
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4.3 Equations of motion 

Rewriting motion equations for the particular case A4 = ttiq = const ., we obtain: 

Proposition 4.1. (Equations of motion: usual particles, four-dimensional observer.) 

Let g be the Minkowski metric, (f> a covariant vector field, and a function on M^. Extremals of 
a Lagrangian system defined by the Lagrangian A = Lds, 



(4.17) L = ~y{u, u) + ^{u) - ^ = -\{{q'f - + i^^- - V'' 

on M X TW^ . and subject to the constraint 
(4.18) 



\ 



1=1 



are curves 7(5) = (s, g'^js )), satisfying the system of mixed second and first order differential 
equations as follows: (3.22) where M = jtiq (equations of motion), and ccf^ — E (equation of the 
constraint, energy equation). 

Proposition 4.2. (Equations of motion: usual particles, three-dimensional observer.) 

Let (f> (3.301 be a one-form (covariant four-vector field), and a function on the Minko wski s pace- 
time Consider on (M x TW^)\Hq adapted coordinates {s,q'' ,t,v\q^) defined by (3.241 and 

(3.251. Extremals of a Lagrangian system defined by the Lagrangian A = Lds, 

,2 s 



(4.19)L = -i(l-y(^T + e(^A.^-F)'^-^, 



2 V 

and subject to the constraint 
moc 



(4.20) 



1 



are curves c(s) — {t{s) , q'' {t{s))) in R*, satisfying the following system of differential equations: 



Too 



1 - 



grad ip 



(equations of motion), together with the following equation of the constraint (mass equation, resp. 
energy equation): 



(4.22) — = TO, resp. 
ds 



Note that the constraint force along solutions of the equations of motion takes the form 
1 V dip 



(4.23) J^c = - 



Hence, along solutions, J-'c = if v = or ip is a constant of the motion. 
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4.4 Variational principle for time-dependent curves 

Taking into account equations of motion ( 4.21^ , we can ask if they are equations for extremals of 
an (unconstrained) variational principle for curves t — > {q^ (t) , q'^ (t) , (t)) in M.^. This means, we 
ask about a Lagrangian 1-form £{t,q'',v'')dt on M x K'^ such that equations (4.211 would identify 
with Euler-Lagrange equations 



(4^24) |§ 



d dC 



dq^ dt dv'- 



= 0. 



The problem can be tackled directly by applying Helmholtz conditions (necessary and sufficient 
conditions for a system of differential equations be variational "as it stands") [7]. However, the 
answer is known: as shown in [12], equations with the left-hand side as in ( 4.21[ ) are variational if 
the force on the right-hand side is a Lorentz-type force. Hence, we can conclude: 



Proposition 4.3. // the scalar field ^ is a constant then equations (4-. 21) are variational as 
equations for extremals of the Lagrangian C{t, q\ v^)dt onM. x M'^, such that 



(4.25) C 



^ + -Av - 
c 



eV. 



Note an interesting connection between vanishing of the induced constraint force and the exis- 
tence of a Lagrangian, namely J-c = implies that dip/dt — 0, i.e. ip — const. 

Remark 4.4. Within the theory of non-holonomic systems it is known that motion equations of a 
Lagrangian system subject to constraints need not be Euler-Lagrange equations of the constrained 
Lagrangian. This is the case also in our situation. Lndeed, the unconstrained Lagrangian 1-form 



(4.26) A = ids = - 



1 - 



{qy - {-Av-eV]—+il:)ds 



q 



gives rise to the constrained Lagrangian 
Ac = Lc'ds = - (^ 

(4.27) 



-moc 



1 



-Av-\-eV-\- 



mJ 



dt 



= - moc 



1 - 



-Av-\- eV + 



mJ 



\\t 



1 



+ TT^T^nC ds 



defined on the constraint Q+. The relevant part of the Lagrangian Ac for considering time- 
parametrised curves is its horizontal part with respect to the projection onto the time axis, i.e. the 
1-form Cdt, with 



(4.28) C 



-rriQC 



1 



+ - Av - eV - 



However, Euler-Lagrange equations of C do not coincide with equations (4-. 21) unless ijj is a 
constant. 



4.5 Dicke force 

Equations of motion in proposition |4.2| can be equivalently expressed in a way that admits to 
compare the obtained forces with a hypothesis due to Dicke on non-Lorentz type interactions [S] . 
Denote exp ^ — and put 

lb 

(4.29) /i = mo = mo e^. 
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Let us express motion equations (4.21 ) in terms of the new scalar potential fi and "mass" mg. We 
get 



and 

(4.31) ^^ = -e-''^, 
^ ' dt dt' 

hence, 



(4.32) e-''^(^2l=) ^ - moc^ ^1 - ^ grad^. 



'1 

Multiplying this equation by we have 



(4.33) ±(^tL=)^e'^J^L-rhoc'\ll-^-^gr^df,. 



1 - ^ 
V ^ 

However, 

(4.34) grad toq = 1^0 grad e'' = niQe^ grad /i = toq grad /z, 
and we finally obtain 



Proposition 4.5. FFif/i /leZp 0/ i/ie mass Junction rho, depending upon a scalar potential ip on . 
and defined by 



(4.36) mo = mpexp 



equations of motion in proposition \4T^ take the following equivalent form: 
• in absence of the electromagnetic field: 



• in presence of the electromagnetic field: 



(4.38) A(^2^) = !!^^,_,^Ji_!^gradmo. 
dt\ _ v^/ mo V 



Hence, the particle moves like having a non-constant rest mass and subject to a force depending 
upon the particle speed and the gradient of the mass. The influence of an electromagnetic field on 
the particle's dynamics depends not only upon the particle charge, but also upon exp ('0/(moC^)) = 
mo /mo. 
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Let us denote 



(4.39) To = \/ 1 - gradmo 



and call J-jj Dicke force. 



Remark 4.6. The obtained formula for the force Td complies with a prediction of a relativistic 
non-Lorentz type interaction due to Dicke [5 . Although in his paper a formula or, at least, exact 
arguments are missing, Dicke conjectured the existence of a force, originating from the mass distri- 
bution of the universe, and acting on particles with nonzero (positive) rest mass. By his hypothesis, 
the particle should move like having a nonconstant mass depending upon a scalar field, and subject 
to a force depending upon the particle 's nonconstant mass, and upon its speed proportionally to the 
factor — -^/l — jf? . Remarkably, this conjecture was influenced by Brans and Dicke' s modifica- 
tion of general relativity, published earlier, attempting to adapt the general relativity theory to the 
Mach principle (see 

5 Tachyons 

Let us suppose < and constant on W^. We set 
(5.1) M = -ml, 

where mo > 0. Particles of this kind are particles with negative square of mass, called tachyons. 

Again we consider either fibred coordinates {s,q'^,q"), 1 < a < 4, or adapted coordinates 
{s,q\t,v\q^), 1 < / < 3, on (M X TR*)\Ho. 

5.1 The evolution space 

The constraint Q C (M x TR'^)\Ho is given by equation 



(5.2) (qy-j^iqn'^-mlc^ resp. (g^)^ f 1 - ^ 



The evolution space is then the submanifold 



2 2 

-mnC . 



(5.3) q^ = 



1 - ^ 



For tachyons one can introduce instantaneous mass and (kinetic) energy by the following for- 
mulas 

/ ■ "^0 *'7io 

(5.4) m = . I.e.: m = — , where v > c, or m = — , where v < c. 



(5.5) S^mc^. 

By the above, tachyons can be regarded either as particles with real positive mass mo moving with 
a speed v > c, or as particles with imaginary mass irriQ, moving with speed lower than the light 
(like "normal particles"). 
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Figure 2: Projection of the evolution space of a tachyon {Ai < 0). 

In our setting, when G M, tachyons appear as particles with real mass mo > moving with 
a speed greater than the light speed. In this case the quantities m and £ are real. Note that if the 
speed of a tachyon is increasing, its mass m is decreasing (t; — > oo means that m —> and £ —> 0, 
and conversely, v ~* c means that to ^ oo and £ oo), and that 

(5.6) TO = Too if V — cV2. 

Similarly as in the case > we can "visualise" the constraint Q and the evolution space 2+ 
in the "velocity space". Supressing the q'^ dimension and denoting again {q^,q^,q'^) = c{x,y,w), 
the manifold Q is given by equation 



(5.7) 



(Fig. 2). 



y 



-1 



TOn 



TOn 



5.2 Four- momentum and energy 

We can introduce four-momentum in a full analogy with the A4 > case setting 
(5.8) p^g{u,-). 

In coordinates p — Padq'^ = piAq^ +p4di, and we obtain the same formulas for the components as 
in the previous section: 

{Pa)l<a<A = ( , Cq ), {p )l<a<A = ( , Cq ). 

~ ^ c ^ ^ c 

For the equation of the constraint Q, i.e. g{u,u) ~ —rn^c^, we obtain 



(5.10) PaP" - -mlc^ 



or pip'- + ^^4^"* = -TOqC^. 



On the evolution space 2+, where 



(5.11) = 



TOqC 



£ 
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components of the four-momentum take the form 

(5.12) (p^)i<o-<4 = (-mw, £/c), {p'')i<„<i^ {mv, E /c). 

Finally, with help of ( 5.10| we can compute the energy: 



(5.13) £ — —pip^ — rn^c? = c^Jp^ — m^c^, 

where, as in the previous sections, — p ■ pis the usual square of length of the three-vector p. 
5.3 Equations of motion 



Rewriting motion equations (3.22) for M ~ ^"t-q = const ., we obtain: 

Proposition 5.1. (Equations of motion: tachyons, four-dimensional observer.) Let g 

he the Minkowski metric, 4> covariant vector field, and Tp a function on M*. Extremals of a 
Lagrangian system defined by the Lagrangian A ~ Lds, 

3 

(5.14) L - -yiii, it) + <j>{u) - ^ = ((g^)' - Y.^q'f) + - ^' 
on M X TM'*, and subject to the constraint 

(5.15) q^= 5](g')2-m2c2 

\ 1=1 

are curves 7(5) = ( s, q'^{s )), satisfying the following system of mixed second and first order dif- 
ferential equations: (3.22) where M. — —ttiq (equations of motion), and cq^ — £ (equation of the 
constraint, energy equation). 

Proposition 5.2. (Equations of motion: tachyons, three-dimensional observer.) Let 4> 

(3.301 be a one-form (covariant four-vector field), and ip a function on the Minkowski spac e-tim e 
M^. Consider on (M x TR^)\Ho adapted coordinates {s,q^ ,t,v'' ,q'^) defined by (3.24) and (3.25). 
Extremals of a Lagrangian system defined by the Lagrangian A = Lds, 

C 



(5.16) L^-\{l~^^)iqr + e{lAv-Vy--i., 



and subject to the constraint 
(5.17) 



are curves c(s) — {t{s) , q\t{s))) in M'*, satisfying the following system of differential equations: 



mov 



mo 



1 grad'0 



(equations of motion), together with the following equation of the constraint (mass equation, resp. 
energy equation): 



(5.19) 



dt 
ds 



m, resp. 



c^^^£. 
ds 
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The constraint force along solutions of the equations of motion in this case takes the form 
1 V 

(5.20) J^c 



moc^ _ ^ dt 

Along solutions, J-c if u = or ■0 is a constant of the motion. 



Proposition 5.3. If ip — const., equations of motion {5.18) are variational as equations for 
extremals of the Lagrangian £{t, q\ v^)dt on M x M? , such that 

(5.21) C = moc^\h^ - 1 + -Av~eV. 

V c 

Finally, let us investigate a Dicke-type influence on tachyons. Denote exp ijl = and put 

(5.22) /i = 2"^, mo = TOoe^. 

An analogous procedure as in the case of "normal particles" directly leads to the following result: 

Proposition 5.4. With help of the mass function mg, depending upon a scalar potential ip on M*, 
and defined by 



(5.23) rriQ — niQexp 



nlc 



equations of motion for tachyons take the following equivalent form: 
• in absence of the electromagnetic field: 



• in presence of the electromagnetic field: 



(5.25) — , — Tl+C' \ —~1 grad ttiq. 

dt V /y2_ niQ V 



We can see that in case of particles with negative square of mass Dicke force takes the form 



(5.26) -^^n = V - 1 grad Too. 



6 Particles with zero mass 



It remains to study the case A4 — 0, i.e. particles with zero mass. Let us again consider on 
(M X TR*)\Ho adapted coordinates (s, g', i, w', (j''), 1 < Z < 3. 
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(6.1) (g4)2_^(^P)2^o, resp. {q^y fl - ^) = 



Figure 3: Projection of the evolution space of a massless particle. 

6.1 The evolution space 

The constraint Q C (M x TM'*)\7io is now given by equation 

p=l X 

Since ^ 0, it holds v = c, i.e., particles of this kind move with the light speed. 

Note that in coordinates {s^q^t^v^ ,q^), the constraint Q+ cannot be expressed in the form 
— h{s,q'' ,t,v''), so that we have to represent the evolution space in the form 

(6.2) Q+: v^c, > 0. 

Drawing the constraint Q and the evolution space Q+ in the "velocity space" similarly as in the 
previous two cases, i.e. supressing the q^ dimension and denoting {q^,q^,q^) = c{x,y,w), we get 
for Q 

(6.3) w'^ -x^ = Q 
(Fig. 3). 

6.2 Four-momentum and energy 

Since our definition of four-momentum is universal for all particles, we have p = g{u, •), and in 
coordinates 

/ \ / -1 -2 •4\ / •4\ 

{Pa)i<a<i = ("q ,~q ,q ) = ( , q ), 

(6.4) " " ^ 

(p<t)i<<t<4 = ( , cq'^), (K)i<<T<4 = ( , eg"). 

^ ^ c ^ ^ c 

For the equation of the constraint Q, i.e. g{u,u) — 0, we obtain 



(6.5) Pap" = 0, or pip^ + -^PiP"^ = 0, i.e. p^p'^ - c^p^ = 0, 
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where, as usual, = p- p = 

On the evolution space it holds v = c. Introducing the unit vector in the direction of v, 

i.e., 

(6.6) V = ccv, 

we obtain the components of the four-momentum on as follows: 

(6.7) {Pa)l<a<4 = q'^{-ev, C), {p" ) l<a<4 ^ {Sy , c) . 

Let us introduce energy by the same formula as in the M ^ cases: 

(6.8) £ = cq^. 

Then the equation of the constraint takes the form £^ — c^p^ — 0, £ > 0, hence, on the 
evolution space, 

(6.9) £^cp. 

6.3 Equations of motion 

It is to be stressed that the case = is much different from the above two, when A4 0. The 
reason is that the constrained euqations of motion are singular, in the sense that they cannot be 
put into the normal form. Recall that they read 

(6.10) - . ' 

904 d(j)i \ 9-0 dip q 



where 1 < ? < 3, for dynamics proceeding in the evolution space Q+. r] 

Let us express these equations in adapted coordinates {s,q\t,v,q^). Since along the con- 
strained paths we have 

'^i' 1 d(z;J g4) q^ ( ^ dq^ .4 dt;^" 

d7 = c^^d^^^r dr + ^ ^ 

we obtain after some calculations 

^^^^^ C c2 j di ""^ \dq^ dql)^ c2 \dt ""dqjJ^Vdq^ dt 

1 f dtp ^ dtp\ 
I ' 

where we have used the constraint conditions v = c and > 0. 



94 ' c dt 



^As it is known, singular equations do not obey the Newton's determinism principle, i.e., by initial conditions in the 
evolution space a unique solution need not be determined. Lagrangian dynamics of this kind are studied by the Dirac's 
theory of constrained systems, or by direct methods developed in [11] (see also [12] and [17]). The latter methods are 
suitable also for systems that are not Lagrangian, as in our case. 
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With notations {411)1=1.2,3 ~ f^, (f>4 — ~f^i ^^^d the unit vector e„ in the direction of v, the 
motion equations take the vector form as follows: 



(6.13) 



~dt 



— e e„ X rot A grad V + \ e„ 

\ c at 



I OA 
~c~dt 



grad V]]ey 



dt 



(above, the dot means the scalar product of two three- vectors) . However, these equations can be 
simplified. First, we notice that Ci, • = 1, hence 

(6.14) iA(e- .e.) = e„.^ = 0. 



Next, denote 



(6.15) T^e ( e; X roti*- - gradT^ 

\ c at 



1 dA 

elk 



grad V 



then obviously 

(6.16) f-v^f-ey = 0. 

Now, applying to the motion equation the scalar product with e^, we immediately obtain 

(6.17) ccy ■ grad-0 + — ^ v ■ gradV" + ~g]^ ^ ^ ^■ 



(Note that the same result provides eq. (3.431 where = 0) 



Finally, we notice that the following fields compatibility condition 
(6.18) J?-gradi/; = 

has to be satisfied. Indeed, due to the above, equations 



de 

(6.19) 9^37" = -^^j and ceu-grad'0 



= 



dt ' " ^ ^ dt 

have to be satisfied simultaneously. Differentiating the latter with respect to t we arrive at 



(6.20) 



P ^ ^ de^ ^ grad ip + ce^-^ grad V' + t; ^ = c ^ • grad ip + ce^- grad *^ + ^ *^ 
dt dt dt at dt dt at dt 



d&„ 
"d^ 



jradi/) 



9" 



■ J- ■ grad V', 



proving our assertion. 

Summarising the results, we can conclude: 



Proposition 6.1. Lef0(3.3Ol he a one-form (covariant Jour-vector field), and ip a function on the 
Minkowski space-time M'*. Consider on (M x TM.'^)\Hq adapted coordinates {s^q^t^v^ ,q^) defined 



by (3.241 and (3.25). Extremals of a Lagrangian system defined by the Lagrangian A = Lds, 

,2 s 



(6.21) L 
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and subject to the constraint 



(6.22) v = c, q^>Q, 

are curves c(s) — {t{s) , q'' {t{s))) in M"*, satisfying the following system of differential equations: 



(6.23) 



f de„ / ^ J- 1 

— = e e„ X rot A — - grad V + I e, 

c at \ c ot 

at ot 



1 dA 
~cl)t 



grad V]]e^ 



together with the following equation of the constraint (energy equation): 

(6.24) c^^=£>Q. 
as 

As a consequence of these equations, the motion is allowed only if the vector fields T and grad "ifi 
are orthogonal (fields compatibility), and proceeds in such a way that the vector fields v and dv/dt, 
T and v, and grad^/i and dv/dt are always orthogonal, where v = ce^, and T is the force on the 
right-hand side of the first equation. 

Let us study the motion equations in more detail. 

Assume the electric charge e be zero. Then T = t) and the motion is governed by equations 



(6.25) 



de;, 



0, 



d^ 
d^ 



C 67, 



^rad i]} 



0. 



The first equation gives us that = k = const ., i.e. the particle moves in along a straight 
line (with the speed of light). More precisely, every straight line in M"^ is an allowed trajectory, 
depending upon initial conditions. Thus, it must hold dip/dt — along every straight line in M.^, 
and consequently, dip/dt — 0, i.e. ^p = const, on M.^. 

However, this result means that the above two equations are incompatible. A release is to 
conclude that massless particles do not feel the scalar field ijj. On the other hand, we have seen 



earlier (cf. e.g. eq. (3.35l) that particles with non-zero mass do feel the field ip. Hence, it turns 



out that every particle possesses a scalar field charge depending upon the particle 's mass, such that 



for massless particles the scalar field charge is zero. With regard to eq. (3.351 we may assume this 
charge being mg, and write 

(6.26) 'ip — niQip. 

Now, proposition |6.2| can be reformulated as follows: 



Proposition 6.2. (Equations of motion: particles with zero mass, 3-dimensional ob- 
server.) Let (p (3.301 be a one-form (covariant four-vector field), and ^ a function on the 
Minkowski space-time . Consider on M x rM*\7Yo adapted coordinates (s, 5', t, v\ g^) defined by 
(3.241 and (3.25). A massless particle is described by a Lagrangian A = Lds, 



(6.27) L 



i('4)('*''+''(;-*"->')7-''*> 

where the charge a — Q, subject to the constraint 
(6.28) v = c, q^>Q. 
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Trajectories are curves c(s) = (t(s),g'(t)) in M^, satisfying the following system of differential 
equations: 



(6.29) 



- — = e le^xrotA--— - grad ^ + (^i- ' ^ + S'^^^ ^ 

c^^^f >0, 
as 



where £ is a function on M*, having the meaning of the kinetic energy of the particle. 

As a consequence of these equations, the motion proceeds in such a way that the vector fields 
V and dv/dt, and T and v are always orthogonal, where v — cCy, and T is the force on the 
right-hand side of the first equation. 



Note that from the form (6.29 1 of the equations of motion it is exphcitly clear that equations for 
massless particles are indeed singular: they do not provide a uniquely determined motion unless 
energy of the particle is fixed. In more detail, 

(i) uncharged massless particles in a scalar field move along straight lines in M'^, with the speed of 
light, and positive energy. 

(ii) Charged massless particles in a scalar and electromagnetic field move with the speed of light, 
the direction of the motion in K'^ is not influenced by the scalar field, is always orthogonal to the 
electromagnetic force, governed by equation 

(6.30) - - = e [e. x rot A - - - - grad + (e. • ( - - + grad v) ) e„ j , 
and depends upon a choice of the particle's kinetic energy £ > 0. 

In both the cases, a trajectory c(s) in can be uniquely determined by initial conditions if energy 
£ (as a function on ]&■*) is chosen. 



A brief summary and conclusions 

We proposed a new approach to foundations of relativistic dynamics. It is based on a treatment 
of a particle as a Lagrangian system subject to a nonholonomic constraint, compatible with the 
special relativity theory. We explored a geometric setting to nonholonomic mechanics, proposed by 
one of us [13] , where the constrained system is considered as "living" on the submanifold defined 
by the constraint: The constraint submanifold thus becomes a genuine evolution space for the 
constrained system. 

The key idea is that relativistic particles of any mass moving in an electromagnetic field (defined 
by a four potential ^) and a scalar field (defined by a scalar potential ?/;) can be described by 

(1) the Lagrange function 

(6.31) L = -\^gav(f(t + e4>„cf - a-^, 

defined on M x TM'', where g is the Minkowski metric on M"', g ~ (—1, —1, —1, 1), and e, resp. a, 
is the particle's charge corresponding to the electromagnetic and scalar field, respectively, 

and 

(2) the nonholonomic constraint 

(6.32) g^.r?" - Mc^ 
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in M X TM"', where Al is a function on the space-time, representing the particle's mass. 

As a main result, we found equations of motion as they appear to a four-dimensional, and to 
a three-dimensional observer. 

We discussed in detail the cases when A4 = const . = mg (usual particles), Ai = const . = — TOq 
(tachyons), and = (massless particles). 

As expected, massless particles move with the light speed, and particles with a positive square 
of mass move with a speed lower than the speed of light; the constant mp has the meaning of rest 
mass. On the other hand, particles with a negative square of mass (tachyons) are particles with 
real positive mass mg, moving with a speed greater than the light speed. The constant toq is then 
the mass of a tachyon moving with the speed c\/2. (Alternatively, tachyons can be understood 
as particles with imaginary rest mass irriQ moving with a speed lower than the speed of light, if 
modeled on a complex Minkowski space). 

It turned out that the scalar field charge a is closely connected with the mass of the particle: 
for massless particles a = 0, and a 7^ in the other cases. This means that the presence of a scalar 
field does not effect the motion of massless particles. On the other hand, all massive particles 
are influenced by a scalar field. We have found that the corresponding force takes the form for 
"normal" particles, respectively for tachyons: 



— * / u'^ — ' / u'^ 

(6.33) J-D = — c Y 1 ^ gradmo, respectively J-d = c y ~^ ~ 1 gradmp, 

where 

(6.34) mo = mpexp ( — ^-^ ) , respectively, mg = moexp ( — 

and we have put a = mg. 

We called !Fd Dicke force in honour of Dicke who predicted existence of a force of this kind [5]. 

The charge e has the usual meaning, however, the corresponding force acting on massive and 
massless particles is different. For massive particles we get the usual formula for the Lorentz force 
(the same for "normal" particles and tachyons). For massless particles the corresponding force 
takes the form 

(1 dA ( 1 dA \ \ 

e; X rot ^ - - — -gradT^ + e; • (^- — +gradV"j j . 

It is important to stress a striking difference between dynamics of massive and massless par- 
ticles. Massive particles are regular mechanical systems in the sense that the dynamics in the 
evolution space obey the Newton's determinism principle (motion equations can be put into the 
normal form, i.e. every trajectory is uniquely determined by initial conditions). On the other 
hand, massless particles are singular mechanical systems: the dynamics in the evolution space do 
not obey the Newton's determinism principle (motion equations cannot be put into the normal 
form, i.e. trajectories are not completely determined by initial conditions). To get a concrete 
trajectory in the space-time M"', one has to fix initial conditions and the energy of the particle. 

We also introduced the four-momentum of a particle as the one form, canonically related with 
the 2-tensor field g. With help of four momentum p = (pa), the evolution space (— the constraint 
submanifold) is given by equation 

(6.36) PaP"^ — J^qC^, respectively, PcrP"^ — —itiqC^, respectively, PaP"^ = 

for "normal" particles, tachyons, and massless particles, respectively. On the evolution space 
components of the four-momentum were shown to be as follows: 

Ev E 

(6.37) {pc) = {-mv, £/c), respectively, (pcr) = ( ^, £/c) = -(-e^, 1) 

c 
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for massive, respectively, massless particles, where 



(b.ijoj m— — , respectively, m- 



for "normal" particles, respectively, for tachyons, and E is the kinetic energy of the particle, defined 
(in all the cases) by formula 

(6.39) E = c(f. 

Consequently, on the evolution space, energy is given by the following formulas: 



(6.40) £ = + rn^c^ , respectively, £ — ~ ttiqC^ , respectively, £ = cp, 

for "normal" particles, tachyons, and massless particles, respectively, where p — \Jp ■ p is the usual 
length of the momentum three- vector p. 

Finally let us mention variationality aspects of the obtained equations of motion. In the 
particular case of "usual" particles with positive square of mass moving in an electromagnetic 
field, the "(3+ l)-dimensional" equations of motion we have obtained, take the usual form. Hence, 
as it is known, they can be obtained from a variational principle as Euler-Lagrange equations of 



the Lagrangian (1.1). On the other hand, rather surprisingly, in presence of a scalar field, the 
obtained equations do not come from an (unconstrained) variational principle^ In particular, it 
might be interesting that trying a "natural" extension of the Lagrangian ( |1.1| by adding to the 



Lagrangian a scalar field term, does not provide the correct equations of motion (recall Sec. 4.4) 
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